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Abstract 

A systematic treatment of the magnetic fluctuations effect on the properties of 
the normal-to-superconducting phase transition in a zero external magnetic field is 
given within the self-consistent approximation and the quasi-macroscopic Ginzburg- 
Landau model. New results for thin superconducting films are presented. Thermo- 
dynamic quantities having a direct experimental interest as the order parameter 
jump, latent heat, and specific heat are considered and numerically evaluated for 
bulk Al and thin Al films. The possibility for an experimental verification of the 
theoretical predictions is discussed. 

1 Introduction 

In 1974 Halperin, Lubensky and Ma (HLM) pP showed that the magnetic fluctuations 
change the order of the superconducting phase transition in a zero external magnetic field 
(Ho = \Hq\ = 0)), i.e., the order of the phase transition from normal-to-uniform (Meissner) 
superconducting state at T c0 = T C (H = 0). In the mean-field approximation, when both 
magnetic and superconducting fluctuations are neglected, this phase transition is of second 
order; see, e.g., Refs. [3 El- Moreover, the fluctuations ip(x) = [ip(x) — (ip(x))] of the su- 
perconducting order parameter ip(x) towards the statistical average (ip(x)) are extremely 
small and can be safely ignored in usual low-temperature (T c0 < 20 K) superconductors. 



For a long time these superconductors have been considered as an excellent example of 
a standard phase transition of second order described by the mean-field approximation. 
When the magnetic fluctuations are taken into account in the Ginzburg-Landau (GL) 
free energy F(ip,A) of superconductor [2j, the same normal-to-superconducting phase 
transition in a zero (mean) external magnetic field (H = 0) is found to be a weakly- 
first order phase transition with a very small latent heat which cannot be observed by 
available experimental techniques pp. The effect of a magnetic fluctuation change of the 
superconducting phase transition order, called HLM effect, is very weak in bulk (three 
dimensional, or 3D) superconductors even in Al where the GL number k is very small 
(k < 1) - a circumstance which is in favor of the effect PIE]. 

In this paper we shall investigate this fluctuation-induced first order phase transition in 
thin (quasi-2D) superconducting films. Bulk superconductors will be also discussed in 
order to compare them with the behavior of thin films. We shall use a self-consistent 
approximation jl], in which the fluctuations 5tp of ip are neglected but the magnetic 
fluctuations are completely taken into account. Note, that the so-called tree approxima- 
tion p| does not yield the HLM effect and the self-consistent, or mean-field-like, approx- 
imation, mentioned above, is the simplest analytical method for an investigation of this 
phenomenon. 

The present study is intended to provide enough theoretical results about the behavior of 
measurable physical quantities directly related to the phase transition properties and in 
this way to ensure a theoretical basis for future experiments on the existence of the HLM 
effect. The need of an experimental observation of the HLM effect is very important be- 
cause the effect remains a theoretical paradigm without a reliable experimental verification 
although its mechanism - the interaction of gauge fields in a quite universal Abelian-Higgs 
model - is of fundamental interest for different fields of physics as pure [3 13 U3 HH U2 UB] , 
and disordered [TU H3 E3 UH UH1 EES 123 121] superconductors, quantum phase transi- 
tions [22122], scalar electrodynamics [21], liquid crystals [231231231211123, and cosmol- 
ogy [S3 EI]- On the other hand, there are some theoretical studies, based on Monte 
Carlo simulations [32], the so-called dual model [33 EH], and certain variants of the 
renormalization-group (RG) J35] I36|. in which no evidence of HLM effect was reported; for 
a discussion of this point, see, the review article [HZ]- Therefore, in the modern theory of 
phase transitions the problem for the existence of HLM effect is controversial and cannot 
be easily solved without a hint from the experiment. The experimental research of the 
effect in liquid crystals cannot be considered reliable although the reported results are in 
favor of its existence. 

Recently, it has been shown [§| that the HLM effect is stronger in quasi-2D supercon- 
ducting films than in bulk superconductors and the preliminary evaluation of the relevant 
physical quantities like the order parameter jump and the latent heat at the equilibrium 
point of the fluctuation- induced first order transition in superconducting films gives for 
them several orders bigger values than for those in bulk materials ^3 E] El EI] • This 
result reopens the problem for an experimental search of HLM effect in type I supercon- 
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ductors, in particular, in thin films of type I superconductors with relatively small GL 
parameter k. Here we shall investigate this problem in a comprehensive way. 

We shall neglect the fluctuations of the superconducting order parameter because their 
effect on the thermodynamics of the superconductor is very weak; see, e.g., Refs. j2HH]- 
Within this approximation, the problems in the scope of the work can be considered 
without the use of RG, as well as of numerous and quite interesting RG results available 
in the literature; for a review, see, e.g., Ref. [3*7] . 

In Sec. 2 we present a derivation of the effective free energy of a D-dimensional supercon- 
ductor. In Sec. 3 we give the first thorough investigation of the effective free energy for 
bulk superconductors. In Sec. 4 the quasi-2D superconducting films and the validity of 
the Landau expansion are discussed. In Sec. 5 we summarize our main conclusions. 



2 EFFECTIVE FREE ENERGY 

2.1 Model considerations 

The GL free energy j2J of a D-dimensional superconductor of volume Vd = {L\...Ld) is 
given in the form 



F(ip,A) = / d u x 



*i 2 + > 4 + ^ 



+ 8^ 



In Eq. (JTJ) the first Landau parameter a = a Q (T — T c0 ) is expressed by the critical tem- 
perature T c0 = T C (H = 0) in a zero external magnetic field (H = \H\), b > is the second 
Landau parameter and e = |e| is the electron charge. The square B 2 of the magnetic 
induction B = (H + AtcM), is given by the vector potential A(x) = {Aj(x), j = 1, D} 
in the form 



B' 



iv 

2 ^ 

i, o = i 



(2) 



here the vector potential A(x) obeys the Coulomb gauge V- A{x) = 0. For a 3D supercon- 
ductor the relation B — [V • A(x)] can be used and when B = Bq is uniform along the z- 
axis, the Landau gauge A (x) = B (—y/2,—x/2,0) can be applied. This representation 
can be generalized for D > 2 - dimensional systems, where the magnetic induction B is 
a second rank tensor: 



Boij = B (5a5j2 — &/2&, 



a, 



(3) 



If we use the notation x = (xi,X2,r), where f is a (D — 2) - dimensional vector perpen- 
dicular to the plane x 2 ), in the 3D case we shall have r — (0, 0, z), and 



(4) 



where e^i is the antisymmetric Levi-Civita symbol. The Landau gauge and Eqs. © - © 
can be used for uniform B = B when SB - fluctuations are neglected; see, e.g. Ref. 0. 
In the prevailing part of our study we shall apply the general Coulomb gauge of the field 
A(x) which does not exclude spatial dependent magnetic fluctuations SB(x). 

In nonmagnetic superconductors where the mean value (M) = (M — 5M) of the magneti- 
zation M is equal to zero in the normal state in zero external magnetic field, the magnetic 
induction in presence of external magnetic field takes the form: 

B = H + 8H(x) + A7r5M(x) , (5) 

where H is the (uniform) regular part of the external magnetic field and 5H is an irregular 
part of H created by uncontrollable effects. We neglect the irregular part 5H and set 
Hq = 0, then B contains only a fluctuation part B = SB(x) = 4nSM(x) that describes 
the diamagnetic variations of M(x) around the zero value (M) = due to fluctuations 
5ip(x) of the ordering field if)(x) above (T > T c0 ) and below (T < T c0 ) the normal-to- 
superconducting transition at T c0 . Note, that the non-fluctuation part A = [A(x)—SA(x)] 
corresponds to the regular part B = (H +(M)) = of B in nonmagnetic superconductors 
((M) = 0) in a zero external magnetic field (H Q = 0). Then we can set A Q (x) = and, 
hence, 5A(x) = A(x), so we have an entirely fluctuation vector potential A(x) which 
interacts with the order parameter ip(x). This interaction can be of type ip 2 A and ifj 2 A 2 
and generates all effects discussed in the paper. 

We accept periodic boundary conditions for the superconductor surface. This means to 
ignore the surface energy including the additional energy due to the penetration of the 
magnetic field in a surface layer of thickness equal to the London penetration depth A(T) = 
Ao|£o| _1//2 , = \T—T c0 \/T c0 ; Ao = (mc 2 b/8Tce 2 aoT c0 ) 1 ^ 2 is the zero-temperature value of A. 
This approximation is adequate for superconductors of thickness L ^> A(T) ^> a , where 
a is the lattice constant and L Q = min{Li, i = 1, ...,£)}. As we suppose the external 
magnetic field to be zero (H = 0) or very small in real experiments, the requirement 
Lq ^> A(T) cannot be satisfied and we take into account only the condition Lq ^> ao- 

In microscopic models of periodic structures the periodic boundary conditions confine the 
wave vectors ki = {ki = (27T7ij/Lj); % = 1, ...,£)} in the first Brillouin zone [— (n/a ) < 
ki < (7r/ao)] and the expansion of their values beyond this zone can be made either by 
neglecting the periodicity of the crystal structure or on the basis of the assumption that big 
wave numbers k = \k\ have a negligible contribution to the calculated quantities. The last 
argument is widely accepted in the phase transitions theory where the long-wavelength 
(kao -C 1) limit can be used. In particular, this argument is valid in the continuum limit 
(Vp/a® — * oo). Therefore, for both crystal and nonperiodic structures we can use a cutoff 
A ~ (vr/ao) and afterwards to extend this cutoff to infinity provided the main contributions 
in the summations over k come from the relatively small wavenumbers (fe < A). Note, 
that here we make a quasimacroscopic description based on the GL functional © which 
means that the microscopic phenomena are excluded from our consideration. 

The GL free energy functional takes into account phenomena with characteristic lengths 



4 



£o and Ao or larger (£ and A) where A(T) is the London penetration length mentioned 
above and £(T) = £oN~ 1//2 is the coherence length 0; here £ = (^ 2 /4maoT c0 ) 1//2 is the 
zero-temperature coherence length. In low-temperature superconductors £o and Ao are 
much bigger than the lattice constant ao- Having in mind this argument we shall assume 
that in our investigation A <C (ii/a ). Whether the upper cutoff A is chosen to be either 
A ~ l/£ or A ~ 1/Ao is a problem that has to be solved by additional arguments (see 
Sec. 3.3). 



We shall use the Fourier expansion 



A 3 {x) 



V 2 

1 D 



and 



lK*) = -VZ>® e **' (7) 

where the Fourier amplitudes Aj(k) obey the relation A*(k) = Aj(—k) and k.A(k) = 0. 
The Fourier amplitude ip(k) is not equal to ip*(—k) because ip(x) is a complex function. 
For the same reason ip(0) = ip(k = 0) is a complex number. 



2.2 Approximations 

The total ignoring of both superconducting and magnetic fluctuations in Eq. (JTJ) leads to 
the familiar free approximation where the GL equations [2j should be solved. Note, that 
the free, or mean-field, approximation is the lowest order theory within the framework of 
the loop expansion, e.g., see [HI EE]- The systematic treatment of the fluctuation effects 
in the asymptotic vicinity of the phase transition point can be given by RG. 

The effect of the superconducting fluctuations 5ip(x) on the phase transition properties is 
restricted in a negligibly small vicinity (|T — T c o| ~ 10~ 12 -j- 10~ 16 ) of the temperature T c q 
and we shall assume that Sip(x) = 0, i.e., ip ~ (ip(x)); from now on we shall denote (ip(x)) 
by ip. So we apply a mean-field approximation with respect to the order parameter ifj(x). 
Within this approximation we shall take into account the 5A(a;)-fluctuations for Bq = 0, 
i.e., A(x) = 5A(x). Furthermore, the A (^-fluctuations can be integrated out from the 
partition function, defined by: 

Z{i>) = J VAe- F ^ A)/kBT , (8) 
where the functional integral J T>A is defined by 

/ao D 
J] J] dA^Sidivlix)} . (9) 
00 3=1 xGV D 

The integration is over all possible configurations of the field A(x)\ the ^-function takes 
into account the Coulomb gauge. 
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The partition function Z(ip) corresponds to an effective free energy T 



<F. D 



k B T\nZ(ij) , 



(10) 



The magnetic fluctuations will be completely taken into account, if only we are able 
to solve exactly the integral (|5|). The exact solution can be done for a uniform order 
parameter if). The uniform value of if) is different from the mean- field value of if) because 
the uniform fluctuations of if){x) always exist, so we should choose one of these two 
possibilities. The problem for this choice arises after the calculation the integral (JEJ) 
at a next stage of consideration when the effective free energy Tv is analyzed and the 
properties of the superconducting phase {if) > 0) are investigated. The effective free 
energy is a particular case of the effective thermodynamic potential in the phase transition 
theory [31 EH] and we must treat the uniform if) in the way prescribed in the field theory 
of phase transitions. It will become obvious from the next discussion that we shall use 
a loop-like expansion which can be exactly summed up to give a logarithmic dependence 



Because of the spontaneous symmetry breaking of the continuous symmetry in the ground 
state, the ordered phase if) > 0, i.e., the effective free energies discussed in this paper 
depend on the modulus \ip\ of the complex number if) = \if>\e ie but not on the phase 
angle 9 which remains arbitrary. That is why we shall consider the modulus \if)\ as an 
"effective order parameter" because the angle 9 does not play any role in the phenomena 
investigated in the paper. The quantity \if)\ remains undetermined up to the stage when 
we define the equilibrium order parameter \ip \ by the equation of state [dF^ip) / dip] = 0. 
This equation gives the equilibrium value if>o of if) and the difference 5ipo = {ipo - if)) 
can be treated as the uniform (zero dimensional) fluctuation of the field if){x). The in- 
dependent fluctuations Sip{x) have been neglected because of the uniformity of if). The 
solution ip will be stable towards the uniform fluctuation 5 i/j provided the same solution 
ip = lipole 100 corresponds to a stable (normal or superconducting) phase; the phase angle 
#o remains unspecified. Therefore, we begin our investigation setting if) uniform but at 
some stage of consideration we shall also ignore the uniform fluctuation Sip and deal 
only with the equilibrium value ipo of if). The equilibrium value will be calculated after 
taking into account magnetic fluctuations, so it will be different from the usual result 
l^o I — (l a l/fr) 1//2 13 when both magnetic and superconducting fluctuations are ignored. 
This simplest approximation for the equilibrium value of if) is obtained from the GL free 
energy provided e = and the gradient term is neglected. Hereafter we shall keep 
the symbol \if>o\ for the equilibrium order parameter in the more general case when the 
magnetic fluctuations are not neglected and shall denote the same quantity for e = by 



The above described approximation neglects the saddle point solutions of GL equations, 
where (ip{x)) is x-dependent. Therefore, the vortex state that is stable in type II su- 
perconductors cannot be achieved. This is consistent with the choice of a zero external 
magnetic field, where the vortex state cannot occur in any type superconductor. These 




V=\Me = 0)\ = {\a\/b) 



1/2 
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arguments can be easily verified with the help of GL equations [2] for a zero external 
magnetic field; the only nonzero solution for ip in this case is given by 77 = (\a\/bY^ 2 
although the magnetic fluctuations A(x) = 5A(x) are properly considered. 

In conclusion we can argue that the described method will be convenient for both type I 
and type II superconductors in a zero external magnetic field, provided the ^-fluctuations 
have a negligibly small effect on phase transition properties T c o = T c (Hq = 0), where T c 
denotes the phase transition line for any H > 0. For type II superconductors in H > 0, 
two lines T cl (H ) and T c2 (H ) should be defined, usually given by H cl (T) and H c2 {T) 0. 



2.3 Derivation of effective free energy 

When the order parameter ip is uniform the functional (JIJ is reduced to 

F(tP,A) = F (ij) + F A (ij) , 

with 



F W = V D (a\ij\ 2 + -\^) 



and 



FaM = ^ / d D x I pWA 2 (x) + \ J2 ( 

I ij'=l ^ 



11 "dAj dA> 



dxi dxj 



(11) 
(12) 

(13) 



Here p = Pol^l 2 an d Po — (8tt e 2 / mc 2 ) . It is convenient to calculate the partition function 
Z(ip) and the effective free energy ^ 7 d('0) in the fc-space, where Eqs. (jHJ) and (fl~3|) take the 
form 



D k<A 



and 
Here 

and 



I f ] J dRzAjfydlmAjfyS k ■ A(k) 
i =1 fc>o 

F A (i/>) = F A (0) + AF A (tfj) . 



AF A (i;) = P J2\Mk) 

note, that we have used the Coulomb gauge k.A(k) = 0. 
Then the partition function (jHJ) will be 

Z{$) = e- FoW/kBT Z A (^) , 

where 

Z A {#) = J VAe- FA ^' kBT 



(14) 

(15) 
(16) 

(17) 



(19) 
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with Fa(?P) given by ()17|) and the functional integration is defined by the rule (JT3J). With 
the help of Eqs. (fTUj) - (fH?j) the effective free energy J r D(i J ) becomes 



where F (ip) is given by Eq. (JT2J) and 

■zty) 



T f {iP) = -k B T\n 



(20) 



(21) 



is the ^-dependent fluctuation part of J-{^). In Eq. ({2D}) the -^-independent fluctuation 
energy {—k B Tln [.2^(0)]} has been omitted. This energy should be ascribed to the normal 
state of the superconductor which, by convention, is set equal to zero. 



Defining the statistical averages 



((•••)) 



jVAe- F ^l k * T {...) 



we can write Eq. (j21|) in the form 

= -k B T\n ( e -^AW/kBT^ 



(22) 



(23) 



Eq. (j2SJ) is a good starting point for the perturbation calculation of !Ff{ijj). We expand 
the exponent in Eq. (|23|) and also take into account the effect of the logarithm on the 
infinite series [3] and obtain in result 

00 f-lV 

^w = E]i(b^ <A ^ W)e ' (24) 

where (...) c denotes connected averages Now we have to calculate averages of the type 

(A a (k 1 ),A /3 (k 2 )...A r (k n )) c . (25) 
Here we shall use the Wick theorem and the correlation function of form 



c\f(k,k') = (M^M-k')) = 5 u ,G?m 



where 



GfAk) = (A t {k)A {-k)) 



Ank B T 



r)- ■ — h h 

^ 2 y '■hj ^iroj 



(26) 
(27) 



and ki = (ki/k). 

The calculation of lowest order terms (I = 1, 2, 3) in Eq. (J24j) with the help of (J25|) - (J27j) 
is straightforward. The infinite series ()24|) can be exactly summed up and the result is 
the following logarithmic function 



1 + 



k 2 



:28) 



The same result for J-"f(i(j) can be obtained by a direct calculation of the Gaussian func- 
tional integral 0. This is done using the integral representation of 5-function in Q 
or (fT4*|) but it introduces an additional functional integration that should be carried out 
after the integration over Aj(x). 

Eqs. ()1L)|). (|2L)jl and (|25|) give the effective free energy density 

f D hl>)=J 7 D^)/V D (29) 

in the form 

/ D W = /„W + A/^), (30) 

where 

/oW0 = aM 2 + ^M 4 (31) 

and 

A/DW= (£z^ Ebl ( 1 + ^. (32) 

Eqs. (12011 and (|29jl - (jH2| are the basis of our further considerations. We should mention 
that the fluctuation contribution Afu^i/j) to f{ip) transforms to a convergent integral in 
the continuum limit 




where Kd = 2 1 ~ D n~ D / 2 /T(D /2) for all spatial dimensionalities D > 2. But the terms 
in the expansion of the logarithm in (}3*2~|) are power-type divergent with the exception of 
several low-order terms in certain dimensionalities D. Therefore, we shall work with a 
finite sum of an infinite series of infinite terms. In our further calculations we shall keep 
the cutoff A finite for all relevant terms in A/d^). This is the condition to obtain correct 
results. 



2.4 Particular dimensions 



For purely 2D superconductor consisting of a single atomic layer, we can use Eqs. (J2HJ)- 
setting D = 2 and calculate Af 2 (ip) with the help of the rule (J3lfy : 



A/2W 



ksT \ \f\2j I/I2N1 ft , A>M 2 \ |,|2i fPo|# 

— j p + p H ) In ^1 + J - PqIVI In (^^- 

The first term of this free energy can be expanded in powers of \ip\ 2 '- 

kp>T\ f , , ( A 2 \ p 2 |V| 4 



(34) 



A/ 2 (V0 



Poi^r+Poi^rin 



Po|# 



2A 2 



(35) 
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Thus we obtain the result from Ref. [123 . This case is of special interest because of the 
logarithmic term in the Landau expansion for f{ip) but it has no practical application for 
the lack of ordering in purely 2D superconductors. 

For quasi-2D superconductors we assume that (2tt/A) > L ^> a , where L is the thick- 
ness of the superconducting film and a more precise choice of the upper cutoff A <C (1/cto) 
for the wave numbers hi is a matter of an additional investigation [Hj (see Sec. 2.1 and 2.5). 
In order to justify this definition of a quasi-2D system we consider the more general case 
of a 3D system of volume V = (LiL 2 L ), where we can take the continuum limit along 
the large dimensions (L\ and L 2 ) of the film because of the assumption L a ^> (27r/A), 
(a = 1,2). The summation over the wave number k = (27ino/L ) cannot be substituted 
with an integration because L <C L a and the dimension L does not obey the conditions, 
valid for L a jJQJ|2JH2j- Therefore, for such 3D system we must sum over fc and integrate 
over two other components {k\ and k 2 ) of the wave vector k. This gives an opportunity 
for a systematic description of the 2D-3D crossover in superconductors [THJ EU 112 EH1 IS] 
which fully justifies the application of more simple treatment for <C L < (27r/A). 

We consider the conditions (2n/A) > L ^> definition of a quasi-2D films of 

thickness Lq. The condition (27r/A) > Lq means that the sum in Eq. ([32]) contains only 
terms with (&o = 0). The summation over k = (ki,k 2 ,0) gives a correct description of 
quasi-2D films of thickness Lq and this can be shown as a limiting case of the more general 
2D-3D crossover described in Refs. [T3*l ITU W2\ |4"3*1 144] . Therefore, for a quasi-2D film we 
have the expression; 

A/(V0 = f A/ 2 (^) , (36) 
where Af 2 (i/;) is given by Eq. (]34|l . 
For the bulk (3D) superconductor we obtain: 
k B T 



A/ 3 (V) 



2tt 



^ ln ( 1 + ^) + ^ M2A _^ fM 3 arctan (_^_ 



(37) 



For the Landau expansion in powers of this form of f^ip) confirms the respective results 
in Refs. 011] and moreover correctly gives a term of type PolV'l 4 which was supposed small 
and neglected in these preceding papers. This problem will be discussed in Sec. 3. 



For 4D-systems A/d(^) becomes 

. / nJW,l 2 \ . . / A 2 \ 1 

(3* 



A%|^|2 + A 4 ln(l + 4^)-^H 4 ln(l ' ^ 



A/4W " 64^ L" rU ' rl A 2 J -U,r, ■ po|# 

The above expression for A/ 4 (■?/>) can be also expanded in powers of \ip\ to show that it 
contains a term of the type |-?/>| 4 ln (y/po\ip\ / A) which produces a first order phase transi- 
tion; this case is considered in the scalar electrodynamics [24 . In our further investigation 
we shall focus our attention on 3D and quasi-2D superconductors. 

The free energy density Af D (i/;) can be expanded in powers of \ip\ but the Landau ex- 
pansion can be done only in an incomplete way for even spatial dimensions. Thus /2(V0> 
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f± (i/j), and f(if>) - the free energy density corresponding to the quasi-2D films, contain 
logarithmic terms which should be kept in their original form in the further treatment 
of the function A/d(V0 m the Landau expansion. We shall do our analysis in two ways: 
with and without Landau expansion of Afr)(ip). These variants of the theory will be 
called "exact" theory (ET) and "Landau" theory (LT), respectively. We shall show that 
these two ways of investigation give the same results in all cases except for quasi-2D films 
with relatively small thicknesses {Lq <C £o)- It seems important to establish the differ- 
ences between two variants of the theory because the HLM effect is very small and any 
incorrectness in the theoretical analysis may be a cause for an incorrect result. By same 
arguments we shall investigate the effect of the factor T in Afj)(ip) on the thermodynam- 
ics of quasi-2D films. This factor can be represented as T = T c0 (l + to) an d one may 
expect that the usual approximation T « T c0 , which is well justified in the Landau theory 
of phase transitions E|, m ay be applied. We shall show for both 3D and quasi-2D 
superconductors, that this way of approximation can be made by neglecting terms in the 
thermodynamic quantities smaller than the leading ones. On the other hand practical 
calculations lead to the conclusion that this approximation cannot be made without a 
preliminary examination because for some quasi-2D films it produces a substantial error 
of about 10%. LT, in which the factor T is substituted by T c0 , will be called a "simplified 
Landau expansion" - SLT. 

2.5 Validity 

The general result (J29|) - (|32j) for the effective free energy f(ip) has the same domain of 
validity [2] as the GL free energy functional in a zero external magnetic field. When we 
neglect a sub-nano interval of temperatures near the phase transition point we can use 
Eq. (|Tj) provided |io| = \T — T c0 |/T c0 < 1, or in the particular case of type I supercon- 
ductors, |io| < ft 2 |2j. Note, that the latter inequality does not appear in the general 
GL approach. It comes as a condition for the consistency of this approach with the 
microscopic BCS theory for type I superconductors 

Taking the continuum limit we have to assume that all dimensions of the body, including 
the thickness L , are much larger than the characteristic lengths £ and A. The exception of 
this rule is when we consider thin films. Especially for thin films of type I superconductors, 
where ((27r/A) > L ^> ao), we should have in mind that £(T) > A(T), so the inequalities 
£ > A > £o > -^o hold true in the domain of validity of the GL theory |£ | < k 2 < 1. 
In Ref. jH] a comprehensive choice of the cutoff A has been made (A = £ ) an d we shall 
discuss this point in Sec. 3 and 4. Note, that the respective conditions for quasi-2D films 
of type II superconductors are much weaker and are reduced to the usual requirements: 
K > 1/V2, \t \ < 1 and (27r/A) > L > a . 

If we do a Landau expansion of /d(V'); m powers of \if}\ 2 the condition /) < A 2 should 
be satisfied. In order to evaluate this condition we substitute \ip\ 2 in p = Pol^l 2 with 
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rf = \a\/b which corresponds to e = (Sec. 2.2). As A 2 (T) = 1/p, the condition for the 
validity of the Landau expansion becomes [AA(T)] 2 ^> 1, i. e., (AA ) 2 3> \t \. Choosing the 
general form of A r = (ttt/^q) where r describes the deviation of A T from Ai = A = (vr/^o), 
we obtain (tttk) 2 ^> |t | ; « = (Ao/£o) is the GL parameter. 

Thus we can conclude that in type II superconductors, where n = (A /£o) > I/v^j the 
condition (p/A 2 ) <C 1 is satisfied very well for values of the cutoff in the interval between 
A = (7r/£o) and A = (7r/Ao), i.e., for 1 < r < (1/k). For type I superconductors, where 
k < l/\/2 the cutoff values A ~ (l/£o) leads to the 'BCS condition (\t \ < k 2 ) for 
the validity of the GL approach. Substantially larger cutoffs (A 3> 7r/£o), for example, 
A ~ (1/Ao) for type I superconductors with k C 1 lead to a contradiction of this BCS 
condition with the requirement /) < A 2 . This inconsistency will be discussed again in 
Sec. 3.3. 

In our calculations we often use another parameter p T = (1/tctk) 2 and, in particular, 
fi = fii = (1/ttk) 2 and in terms of /i the condition for the validity of expansion of 
becomes p\t \ <C 1, or, more generally, p T \t \ <C 1. Choosing r = 1/ir we obtain the 
BCS criterion for the validity of the GL free energy of type I superconductors 0. The 
choice r = (^o/hXq) corresponds to the cutoff A T = 1/Ao- As we shall see in Sec. 3 and 4 
the thermodynamics near the phase transition point has no substantial dependence on 
the value of the cutoff A r but it should be chosen in a way that is consistent with the 
mean-field- like approximation. 

Alternatively, the inequality (p/A 2 ) 1 may be investigated with the help of the reduced 
order parameter <p defined by <p = \ip\frjo, where rj = r](T = 0) = (aoT c0 /6) 1//2 is the so- 
called zero-temperature value of the order parameter within the GL free energy foi^), 
given by Eq. (|3*T|); see also Sec. 2.2. The reduced order parameter ip will be equal to \t \ 
for to < 0, if only the magnetic fluctuations are ignored, i.e., when = rj. Using the 
notation ip, we obtain the condition (p/A 2 ) <C 1 in the form p T <p 2 1. This condition 
seems to be more precise because it takes into account the effect of magnetic fluctuations 
on the order parameter ip. 

3 BULK SUPERCONDUCTORS 

3.1 Free energy 

The effective free energy fz(ip) of bulk (3D-) superconductors is given by Eqs. - (}3*Tj) 
and (J37)) . The analytical treatment of this free energy can be done by Landau expansion 
in small (yfpo~\ip\/ A). Up to order |-?/>| 6 we obtain 

/ 3 (V) « « 3 |^| 2 + - <? 3 H 3 + |h/f , (39) 
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and 



where 

k B TAp 

o 3 = o + , (40) 

Ztv z 

fc = »+^£. (41) 
« = ^ , (42) 

D7T 

The cutoff A in Eqs. (|4*Uj) - is not specified and can be written in the form A T = (vrr/^o) 
as suggested in Sec. 2.5. 

We shall just outline the analysis of the above free energy. It can be shown by both 
analytical and numerical calculations [TD] that |-?/;| 6 -term has no substantial effect on the 
thermodynamics, described by the free energy (|39p. That is why we ignore this term and 
do the analysis in the standard way 3J. The possible phases \ipo\ are found as a solution 
of the equation of state: 

[df(1>)/d\1>\ = . (44) 

There always exists a normal phase \ipo\ = which gives a minimum of fz{ip) for 03 > 0. 
The possible superconducting phases are given by 




\M± = ^ 1±, 1-^ >0- (45) 



160363 
9<7 3 

Having in mind the existence and stability conditions of j^oli-phases |3 , we obtain that 
the |-0o|+-phase exists for (160363) < 9gf and this region of existence always corresponds 
to a minimum of f^{iji)- The |^o| --phase exists for < 03 < (9gf/16&3) and this region of 
existence always corresponds to a maximum of faty), i.e., this phase is absolutely unstable. 
For 03 = 0, \if>o\- = and hence, coincides with the normal phase. For 9gf = (160363) 
we have |^o|+ = |^o|- = (3g 3 /46 3 ) and /3(|^o|+ = h\i>o\-) = (27g|/5126|). Furthermore 
/ad^ol-) > for all allowed values of \ipo\- > 0, whereas 

/ 3 (|Vo|+)<0 for o 3 <(g 3 2 /26 3 ), 

and 

9o 2 

h(\AW) > for (g 3 2 /26 3 ) < o 3 < ' 



I663 

The equilibrium temperature T eq of the first order phase transition is defined by the 
equation f(\ip \ + ) = which gives the following result: 

26 3 (T eq )o 3 (T cq ) = g 2 (T eq ) . (46) 

These results are confirmed by numerical calculations of the effective free energy (|39|) [10 ; 
there also the influence of the |?/^| 6 -term is evaluated. 
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3.2 Entropy and specific heat capacity 



The equilibrium entropy jump is AS = VAs and As = —(df 3 (\ip\)/dT) can be calculated 
with the help of Eq. (|39p and the equation of state 



A S = -|Vo| 2 <Wo|), (47) 
where $(|?/>o|) is the following function: 

The specific heat capacity per unit volume AC = T(dAs/dT) is obtained from (|47j) 

The quantities As(T) and AC(T) can be evaluated at the equilibrium phase transition 
point T eq which is found from Eq. (j4fij) : 



p 3 /2 k B /67i 



T c0 2vr 2 a 6 + (p^ b /2tt 2 A) T c0 V «o 

provided |AT C | = |T c0 — T eq | <C T c0 . Further we shall see that the condition |AT C | <C T c0 
is valid in real substances. The second term in r.h.s. of Eq. (J5U|) is a typical negative 
fluctuation contribution whereas the positive third term in r.h.s. of the same equality is 
typical for first-order transitions [3]. 

To obtain the jumps As and AC at T cq we have to put the solution |V>o|+ found from 
Eq. (g3J) in Eqs. (|47 |) - The result will be: 

3/2 x 2 



A ^-wA a " + ^-[— ) wJ' (51) 

and 



AC= »rK- , -f-l' (52) 



4 ^ Ur„ - f- 4 

where &3 C and g3 C are the parameters 63 and q 3 at T = T eq . As |AT C | = |T c0 — T cq | <C T c0 
we can set T eq « T c0 in r.h.s. of Eqs. (joTj) and (|52"j) and obtain g 3c = q 3 (T = T eq ) « g3(T c0 ) 
and 6 3c « b 3 (T c0 ). 

The latent heat Q = — VT eq As of the first order phase transition at T eq can be calculated 
from Eq. (f5*Tj) . If we neglect the charge (e = 0) which means to set po = <?3 = and 
6 C 3 = b in Eqs. (15 lj) - (|52|) we shall get the result from Ref. PQ for the ratio 

(AT), = ^_ . (53) 
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Here we should mention that Eq. (|52|) gives the jump AC at the equilibrium phase tran- 
sition point of the first order phase transition, described by \ip\ 3 term [3], while AC 
calculated in Ref. [I] is equal to the specific heat jump at the standard second order tran- 
sition AC = (alT c o/b) and is four times smaller. Therefore, we obtain (AT) eq four times 
smaller than the respective value in Ref. p. 

3.3 Numerical values for Al 

In order to do the numerical estimates we represent the Landau parameters ao and b with 
the help of the zero-temperature coherence length £ and the zero-temperature critical 
magnetic field H c0 . The connection between them is given by formulae of the standard 
GL theory of superconductivity 0: £q = / ^ ma oT c o) and if 2 = (inalT^/b). The 
expression for the zero-temperature penetration depth Ao = (hc/2^eH c0 ^ ) is obtained 
from the above relation and Ao = (6/ao^coPo) 1//2 - We shall use the following experimental 
values of T c0 , H c0 and £ f° r Al: T c0 = 1.19K, H c0 = 990e, £ = 1.6/xm, k = 0.01 
The experimental values for T c0 , H c0 and £ vary about 10-15% depending on the 
method of measurement and the geometry of the samples (bulk material or films) but 
such deviations do not affect the results of our numerical investigations. 

The evaluation of the parameters 03 and 63 for Al gives: 

% = (a T c0 ) [t + 0.972 x 10- 4 (1 + t )r] , (54) 

and 

Setting r = 1 corresponds to the cutoff A x = (7r/£ ) (Sec 2.5). For r = (1/k)ai — 10 2 
which corresponds to the much higher cutoff A = (7r/Ao) we have 63 ~ b. i.e., the p\ -term 
in 63, given by Eq. (|41jl. can be neglected. However, as we see from Eq. (|55jl . for r = 1 
the same pg-correction i n the parameter 63 is of order 0.16 and cannot be automatically 
ignored in all calculations, in contrast to the supposition in Refs. [HE]. However, the 
more important fluctuation contribution in 3D superconductors comes from the r— term 
in Eq. ()54|) for the parameter 03. This term is of order 10~ 4 for r ~ 1 and this is consistent 
with the condition |to| < k 2 ~ 10~ 4 but for r ~ 10 2 , i.e., for A ~ (vr/Ao) ~ 10 6 /im, the 
same r— term is of order 10 2 which exceeds the temperature interval (T c0 ± 10 -4 ) for the 
validity of BCS condition of Al (Sec. 2.5). 

These results demonstrate that for our theory to be consistent, we must choose the cutoff 
A r = (vrr/^o), where r is not a large number (r — > 1 -j- 10). To be more concrete we set 
A = Ai = (7r/£ ) as suggested in Ref. [§]. 

The temperature shift t eq = to(T eq ) for bulk Al can be estimated with the help of Eq. (|5T)|) . 
We obtain that this shift is negative and very small: t eq ~ — 10~ 4 . Note, that the second 
term in the r.h.s. of Eq. ()50|) is of order 10~ 4 provided A ~ (l/£o) whereas the third term 
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in the r.h.s. of the same equality is of order 10~ 5 . Once again the change of the cutoff A 
to values much higher than (tt/£o) will take the system outside the temperature interval 
where the BCS condition for Al is valid. Let us note, that in Ref. jTD] the parameter t 
corresponds to our present notation to- But the numerical calculation of the free energy 
function f^ip) in Ref. ^U] was made for the SLT variant of the theory and the shifted 
parameter (t + 0.972 x KT 4 ) was incorrectly identified with t and this lead to the wrong 
conclusion for its positiveness at the equilibrium phase transition point T eq . As a matter 
of fact, the shifted parameter (t Q + 0.972 x 10~ 4 ) is positive at T eq but t eq = t (T eq ) is 
negative. 

Having in mind these remarks, when we evaluate As and AC for bulk Al we can use 
simplified versions of (J51|) and (J52)) which means to consider only the first terms in the 
r.h.s and to take q 3c rs q 3 and b 3c ~ b at T c0 . In this way we obtain 



Q = -T c0 As = 0.8 x 10" 

and 

AC = 2.62 x 10 3 



erg 



K . cm 3 
erg 



cm 3 J 



(56) 



(57) 



The results are consistent with an evaluation of AC for Al as a jump (AC = a^Tco/h) 
at the second order superconducting transition point pQ that, as we mentioned above, is 
four times smaller than the jump AC given by Eq. (J57I) . 

A complete numerical evaluation of the function f 3 {4>) and the jump of the order param- 
eter at T eq for bulk Al was presented for the first time in Ref. ^Hj. The results there 
confirm that the order parameter jump and Q for bulk type I superconductors are very 
small and can hardly be observed in experiments. 

We shall finish the presentation of bulk Al with a discussion of the ratio ([53)1 . It can be 
also written in the form 

<^-?(a)£)'- 

and it differs by a factor 1/4 from the respective result in Ref. pQ. This difference is due 
to the fact that we take AC as the jump at the first order transition temperature T eq 
while in the above cited paper pQ the authors define AC as a hypothetic jump (AC) at 
the standard second order phase transition point. From Eq. ()56|) we obtain 

(AT) C q = 6.7 x 10- 12 (T C 3 ^6 ); (59) 

and multiplying the number coefficient in the above expression by 4 we can obtain Eq. (10) 
from Ref. pQ. 
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4 QUASI-2D FILMS 



Thin quasi-2D films (ao «C Lq < 2tt/A) can be investigated with the help of the respective 
free energy density f{ip) given by Eqs. and (JSTJ, Af D (ip) is taken from Eqs. (jEU) 
and (|36j) . The free energy of quasi-2D superconducting films was derived and analyzed 
for the first time in Ref. """•] using the Landau expansion of A^VO i n powers of \ip\ 2 ; 
see Eq. ("*"""")) . As is shown for the first time by Lovesey jSHj m the simple 2D case the 
fluctuation contribution A / D (■?/;), of form given by Eq. (""""""j). leads to a fluctuation-induced 
first order phase transition. In contrast to 3D superconductors where the first order of 
the phase transition is generated by |?/>| 3 -terni in A/3 (ip), in 2D superconductors the first 
order of the phase transition is a result of the presence of |^| 2 ln|-^| in Eq. ()35|) . But the 
Meissner phase cannot occur in 2D (single atomic layer) superconductors because of the 
strong fluctuations and hence this case is of no interest. In quasi-2D films, where the 
Meissner phase does occur for properly chosen thickness of the film (Lo -C 27r/A) [*""|, the 
change of the order of normal-to-superconducting phase transition is better pronounced 
than in bulk superconductors. This is well illustrated in the above cited paper by 
numerical data for Al films with thickness Lq = 0.1/im. Following Refs. C"T| and the 
arguments presented in Sec. 3.3 we shall choose the cutoff A = tt/^ . 

The expansion of the respective free energy in powers of leads to somewhat clumsy 
analysis and for this reason we shall use the approach in Ref. |TT| IT2"] where the quasi-2D 
films have been investigated with the help of the general form of Af(ip) described by 
Eq. (J211)- In both variants of the theory (ET and LT; see Sec. 2.4) the thickness L of the 
quasi-2D film has an effect on the thermodynamic behavior, that is similar to the influence 
of material parameters cto and b. This is very well seen in the Landau expansion (|3"5|) of 
the free energy f(ip) given by ("ST"*), where the parameters a and b acquire a fluctuation 
contribution that depends on L Q . The influence of L on the thermodynamic properties 
can be considered as a characteristic feature of quasi-2D systems [I2J ESI, a feature, 
absent in purely 2D films jSHl- R is unambiguously demonstrated by several theoretical 
studies of the 2D-3D crossover in systems with slab geometry jT3"l 11*2*1 ESI E3 that the 
L -dependence as given in Eq. ("*""""]) correctly describes quasi-2D films. 

Following Refs. "TT) IT2*] and having in mind the above discussion we can present the free 
energy density f(ip) = (F(ip)/ L1L2) in the form 




(60) 



where 



T iy) = (1+y) In (1 + 2/) hi 



(61) 




(62) 
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In Eqs. (|6(J|) - (|62|) we have set A = (7r/£ ) and introduced the notation = l^l/rfa; 
the quantity i] is defined in Sec. 2.5. Some of the properties of free energy (pjj) were 
analyzed in Ref. JI] for Al films and in Ref. ^2] for films of Tungsten (W), Indium (In), 
and Aluminium (Al). Here we shall summarize and justify the preceding results and, 
moreover, we shall present new results about the properties of the Landau expansion of 
effective free energy. Note, that the function T(y) cannot be fully expanded in powers of 
y because of the term of type (y lny) in Eq. (jrTT]). 

The extensive investigations fTJ^] of films of W, Al and In with thicknesses from 0.05 
fim to 2 /im confirm the intuitive notion that the HLM effect is stronger for smaller 
values of Lq. The numerical analysis shows that type I superconductors with relatively 
small GL parameter k and relatively high critical field H c0 may be the best candidates for 
the experimental observation of the effect. The best material from the above enumerated 
substances seems to be Al; tungsten has an extremely small GL parameter but also a small 
critical field that makes it inconvenient for experiments. The relatively high H c q of In 
results in relatively large latent heat, Q ~ 4.0 (erg/cm 3 ) but for films with Lq ~ 0.05 /xm 
the order parameter jump |^| eq = f e qVo at T eq is twice smaller than that for the respective 
Al films |12j: |V>| eq = 0.05 x 10~ n for In and |^| eq = 0.1 x 10" 11 for Al. We have to stress 
the role of critical magnetic field H c q, a fact established for the first time in [T2| and the 
present paper. 

With the help of data from Refs. (HICTJIEJE] we compared the thermodynamic quantities 
near the first order phase transition point in bulk Al and Al films of L ~ 0.1/iin. They 
are given in Table 1, where t eq = to(T eq ), ip eq = ip(T eq ) is the equilibrium jump of the 
reduced order parameter and \ip\ e q = PeqVo is the order parameter jump at T eq . 

Table 1. Numerical data for bulk Al and Al film of thickness Lq = 0.1//m 



quantity 


teq 


</?eq 


1 


Vieq 


Q (er : 


g/cm 3 ) 


bulk Al 


-0.492 x 10" 4 


0.0032 


Oi 


] x 10 9 


0.8 x 


10~ 2 


Lq = O.lfim 


-0.00147 


0.032 


0.8 


x 10 10 


0.8 
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Figure 1: The function f(<p) for Al films of thickness L = 0.4/iin: the solid line corre- 
sponds to ET, the line of crosses (+) represents LT, and the line of circles (o) stands for 
SLT. All curves are calculated for T eq corresponding to ET (see the text). 

The shift i e q of the equilibrium transition temperature due to magnetic fluctuations is very small 
in both bulk Al and thin Al films so the difference (T eq — T co ) can be neglected in all calculations 
of thermodynamic quantities near T eq . The equilibrium jump ip eCL , or, equivalently, \ip\eq, is one 
order of magnitude higher in the film with Lq = 0.1 /xm than in bulk Al but the latent heat 
Q is 10 2 times bigger for films. These values are almost one order of magnitude higher for 
Lq ~ 0.05 Lira than for Lq = 0.1 Lim. The numerical data in Table 1 are obtained by SLT for 
the bulk Al samples and by ET for the Al film of thickness Lq = 0.1 /im; for the abbreviations 
SLT, LT and ET, see Sec. 2.4. The difference in the numerical results obtained from ET, LT 
and SLT will be discussed in the remainder of the paper. 

The investigation of bulk superconductors yields the same results irrespective of whether we 
analyze the free energy fs(ip) by ET, LT or SLT. The situation in quasi-2D superconductors 
is however different; the three different variants of treatment of the free energy give different 
results, in particular, for relatively small thicknesses (Lq <C £o)- This feature of free energy, 
Eq. is illustrated in Fig. 1, where three curves for three different variants of f(ip) are 

shown for Al films of thickness Lq = 0.4 fj,m. The solid line corresponds to ET, the line 
of crosses represents the LT result, and the line of circles stands for SLT. All three curves 
are calculated for i eq = io(T eq ) = —0.00057, which is the ET equilibrium phase transition 
temperature T eq = 0.9994T c o. Note, that </3 eq is the nonzero global minimum of f(<p) and the 
function f(<p) depicted in Fig. 1 has only one minimum for ip > because all curves are calculated 
in the thermodynamic regime corresponding to the stable Meissner phase. 

The main conclusion that can be made from Fig. 1 is that the two variants of the Landau 
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order parameter 



Figure 2: The function f(<p) for Al films of thickness L = O.lyum: the solid line (-) 
corresponds to ET, the line of crosses (+) represents LT, and the line of circles (o) stands 
for SLT. All curves are calculated for T eq corresponding to the respective variant of the 
theory (see the text). 

expansion give approximately the same quantitative results and therefore, the factor (1 + to) 
in (I6U|) can always be substituted by unity, though the present investigation is intended to quite 
small physical effects. This conclusion is consistent with the argument that allows to use 
the same approximation (1 + to) ~ 1 for the calculation of Q and AC in both variants of the 
theory: ET and LT. Besides, the Fig. 1 shows that both variants of LT give slightly higher 
equilibrium phase transition temperatures T eq and substantially higher equilibrium jumps ip eq 
than ET. Thus, using LT for film thicknesses ao <C Lq < 0.1/iin one may obtain up to 10 times 
higher value of ip eq and up to 10 2 times bigger latent heat Q than the respective values in Table 
1. The problem is whether these higher values predicted by LT are reliable. 

Fig. 2 shows the free energy drawn in the three variants (ET, LT and SLT) of f(<p) for Al films of 
thickness Lq = 0.1/im. In Fig. 2, the curves /(<£>) are drawn at their respective equilibrium phase 
transition points. The variation in t eq (—0. 00148 for the solid line, —0. 00115 for "+" -line and 
— 0. 00046 for o-line) are of order of the typical values of t eq itself so the differences in t eq due to 
the way of calculation cannot be neglected. Although for both variants of the Landau expansion 
(LT and SLT), the quantity ip eq is again practically the same, the difference in t eq is more 
pronounced for smaller thickness of Al film and moreover, both variants of Landau expansion 
are not so good approximation to the result of the exact calculation (ET) as for Lq = 0.4/zm. The 
conclusions, we have already drawn from the results shown in Fig. 1, are completely confirmed 
by the form of the curves from Fig. 2 and, moreover, we see that the deviation of the results of 
LT from those of ET becomes bigger with the decrease of film thickness Lq. 
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Figure 3: The function f((p) calculated from the ET for Al films of thickness L = 0.1//m 
and different cutoffs A: the solid line (-) corresponds to A = tc/Xq, the dashed line 
represents the case A = 1/Ao, and the line of circles (o) stands for A = 7r/£ - All curves 
are calculated at the respective T eq (see the text). 



On the basis of the above observations we may conclude that in all cases when ET and the 
respective Landau expansions give different results, the Landau expansion yields a better estab- 
lished first order transition, with a higher jump </? eq , and hence, bigger values of Q and AC. In 
order to establish where LT is a good approximation we have made systematic numerical calcu- 
lations for Al films of different thicknesses Lq = 0.05 -f 3/um. When Lq is lowered beginning with 
3/um the quantitative differences between the two variants of theory, with and without Landau 
expansion, respectively, become substantial about Lq ~ 0.4/ma. Bearing in mind the condition 
for the validity of the Landau expansion (see, Sec. 2.5) and the requirement for the equilibrium 
jump of the order parameter, < 1 we may suppose that the predictions done with the help 
of the Landau expansions do not satisfy this inequality for Al films with thicknesses Lq < 0.1/zm. 
For this relatively small Lo-size, ET is absolutely reliable. The numerical data show that films 
with Lq > O.lfim are described well quantitatively by the Landau expansion. The differences 
between the curves in Figs. 1 and 2 can be neglected in numerical calculations intended to give 
theoretical predictions for experiments. 

This general result is supported by the following simple argument. In the Landau expansion 
of free energy 1)60(1 for quasi-2D films the parameter a acquires a cutoff (A—) independent 
contribution of the form 

k B Tp 

Aa = -^Lq- ■ (63) 
When we compare Aa with the bare parameter value \a\ = (aoT c o|to|) f° r Al, it is easily obtained 
that | to | will not exceed 10 -4 for thicknesses Lq which are of order l^m or larger. Therefore, for 
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Lq ~ 0.1 jJLva. the Landau expansion gives results which are quantitatively different from those 
obtained by ET. 

We have studied the dependence of the free energy density f(<p) of Al films with Lq = O.l^m on 
the cutoff A. Fig. 3 shows the free energy density /(</?) for three values of the cutoff: A = tt/\q, 
A = 1/Ao, and A = tt/^q. As the cutoff increases from (vr/^o) to tt/Xq ~ 10 2 (7r/£o), the 
equilibrium jump ip eq increases, too. We have already mentioned in Sec. 3.3 that the increase 
of the cutoff A for type I superconductors up to the value (vr/Ao) is inconsistent in the present 
theory. The numerical result for Al films shown in Fig. 3 is, therefore, a demonstration of the 
validity of our arguments about the choice of the cutoff A presented in Sec. 3.3. If we take the 
cutoff A 3> (vr/^o) we shall go beyond the scope of validity of our theory. 

There is a similarity between the breakdown of the present theory for cutoffs A 3> (vr/^o) and 
the breakdown of the condition (p/A 2 ) <C 1 for the validity of LT at small thicknesses Lq. 
In both cases, when there is an inconsistency of the theory, we obtain enhanced values of the 
characteristic jumps of thermodynamic quantities at the equilibrium point of the first-order 
phase transition. 

5 CONCLUSION 

We did a detailed analysis of the HLM effect in bulk (3D-) superconductors and quasi-2D 
superconducting films within the self-consistent approximation introduced in Refs. [Tlllll24j. We 
have studied for a first time the validity of this approximation and calculated thermodynamic 
quantities of direct experimental interest like the equilibrium jumps of the order parameter, 
entropy and specific heat at the point of the fluctuation-induced first-order phase transition 
to superconducting state in a zero external magnetic field. Our investigation is supported by 
numerical calculations for bulk Al and Al films. 

We have presented for a first time a comprehensive analysis of the effective free energy of the 
superconductor in a zero external magnetic field and on the basis of this analysis we compared 
the results from the investigation of the effective free energy without a partial Landau expansion 
with those from the Landau expansion of the effective free energy of quasi-2D films. For quasi- 
2D films of type I superconductors, the Landau expansion leads to reliable results, provided the 
film thickness is above some value depending on the characteristic lengths £o an d Ao, i.e., on 
the material parameters. For Al films the results from the Landau expansion become unreliable 
below film thicknesses Lq ~ 0.1/xm. For quasi-2D films of type II superconductors the Landau 
expansion of the effective free energy can be used for any thickness above the level of destruction 
of superconductivity (Lq ~ 10~ 3 /im). 

Our investigation provides a reliable theoretical basis for a future experimental search of the 
HLM effect in thin films of type I superconductors, where the effect is much stronger than 
in bulk materials. In accord with preceding works [01 110| lll| 1121 113j we have justified earlier 
results which indicate that the HLM effect will be better pronounced in films of materials with 
relatively high values of the critical magnetic field H c q and relatively small thicknesses Lq. We 
cannot be certain which superconducting material provides the best experimental conditions for 
transport or caloric measurements of the jumps of the thermodynamic quantities at the point 
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of the fluctuation-induced first order transition, but from the data for Al, W and In available 
from recent studies [111 I12j and the present paper, the most suitable substance seems to be Al. 
But investigations of other superconductors may put forward materials which are even better 
candidates for the experimental test of the HLM effect. 

Looking for the most convenient material for the experimental search of the HLM effect we 
should have in mind a number of other experimental requirements which are not related to the 
results from the present theoretical investigation. Here we shall briefly discuss the problem of 
the external magnetic field and the possible change of the superconductivity from type I to 
type II with the decrease of the film thickness [Hi]. This change for film thicknesses Lq, which 
are convenient for the experimental study of the HLM effect, is not a great problem because 
we have shown in our investigation that the effect could be observed also in films of type II 
superconductors, provided the external magnetic field Hq is very low so the effect of the vortex 
phase and the magnetic energy jump (Hq/8tt) at the phase transition point is negligible. The 
magnetic energy jump (Hq/8tt) may obscure the HLM effect on the latent heat also in type 
I superconductors and, therefore the experimental problem for the elimination of the residual 
laboratory external magnetic field H$ is common for both type I and type II superconducting 
films. If we take as a basis the latent heat of order 1 [erg/cm 3 ] in Al films with Lq ~ 0.1//m, as 
reported in Sec. 4, the magnetic field which ensures the ratio (Hq/8itQ) <C 1 will be obviously 
about 1 Oe. In thinner films of convenient materials this experimental condition may become 
Ho ~ 10 Oe but no more. 
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